Kinetic Equation for Gluons in the Background Gauge of QCD 
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We derive the quantum kinetic equation for a pure gluon plasma, applying the background field 
and closed-time-path method. The derivation is more general and transparent than earlier works. 
A term in the equation is found which, as in the classical case, corresponds to the color charge 
precession for partons moving in the gauge field. 
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Heavy ion collisions at ultra-relativistic energies are 
widely expected to be a laboratory to study the forma- 
tion and properties of highly excited QCD matter in the 
deconfined Quark-Gluon Plasma (QGP) phase ||]. The 
QGP is considered as a partonic system being at (or close 
to) local thermal equilibrium. Thus, to study the con- 
ditions for the possible formation of QGP in heavy ion 
collisions one needs to address the question of thermal- 
ization of the initially produced partonic medium Q . In 
the theoretical aspect this requires the formulation of the 
kinetic equations involving color degrees of freedom 
and the non-Abelian structure of QCD dynamics. Vari- 
ous models for the initial conditions in ultra-relativistic 
heavy ion collisions suggest that at the early stage the 
medium is dominated by gluon degrees of freedom [Q. 
The kinetic equation for a pure gluon plasma is thus of 
special interest. 

The usual treatment of the gluon transport equation is 
based on the decomposition of the gluon field into a mean 
field and a quantum fluctuation. Under this approxima- 
tion the gluon transport equation then describes the ki- 
netics of the quanta in the classical mean field (|, ^, . 
This picture is somewhat similar to the one used in study- 
ing the energy loss of a fast parton moving in the soft 
mean field [gf To include the classical chromofield 
into QCD in a proper way, one uses the background field 
method of QCD (BG-QCD) introduced by DeWitt and 
't Hooft 0, |ll|, The advantage of BG-QCD is that 
it is formulated in an explicit gauge invariant manner. 

One of the first attempts to derive the gluon transport 
equation in BG-QCD was presented in Ref. But the 
obtained equation is not transparent. The most recent 
work has been done by Blaizot and lancu jl^, ^ in the 
close-time-path (CTP) formalism. There, however, the 
authors focus on formulating the transport equation in 
the vicinity of equilibrium. 
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In this paper, we use the CTP and BG-QCD method 
to derive the kinetic equation for the gluon plasma. Our 
derivation is going beyond previous results as it is quite 
general and transparent. One of the most important fea- 
tures of our work is that a term in the obtained equation 
is shown to correspond to the color charge precession 
term in the classical equation j|, 16, 17[ ^ p^ . 

In the following we use g^^i, — diag(l, — 1, — 1, — 1) as 
the metric tensor, and for elegance of the formula we al- 
ways write Lorentz indices as subscripts and color ones 
as superscripts for the relevant quantities. For the gauge 
field and its strength tensor we denote = A'^T^ and 
Ff,^[A] = F^^^iAlT", where (r°)*J = if^ are the genera- 
tors of the SU{3)c adjoint representation. The two-point 
Green function (GF) or self-energy (SE) are treated as 
matrices, so their color and/or Lorentz indices are some- 
times omitted. 

Applying the background field method, we decompose 
the conventional gluon field into the sum of a classical 
background part A and a quantum fluctuation Q. In- 
cluding the appropriate gauge fixing and ghost terms for 
the background gauge Di^[A]Qj, = 0, the BG-QCD La- 
grangian reads 
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where D''iil[A{x)] — [dxfj. — igAf^{x)Y^ is the covariant 

derivative, C^/C are the ghost/anti-ghost field and a 
is the gauge fixing parameter. 

The above Lagrangian is invariant under the local 
gauge transformation of type I (type II transformation 
is irrelevant to the current problem) where the back- 
ground field transforms as a conventional gauge fields, 
A'^ = UA^U~^+ig~^Ud^U~^, while the gluon and ghost 
fields transform like a matter field, Q'^ — UQ^U~^ W^ - 
Here U{x) = ex.]i[iguj°'{x)T°'] is the transformation ma- 
trix. 

The non-equilibrium dynamics is usually described in 
the CTP formalism |Q. Here the action can be written 
as S'cTP = S+ — S- + K{A±,Q±) where all fields in S± 
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are defined on the positive/negative time branches and 
K(A±,Q±) is the kernel incorporating initial state cor- 
relations. The role of the non-local kernel K{A±,Q±) 
in ^cTP is discussed in Ref. and its coonection 
with the pinch singularity will be discussed else- 
where. The GF for the gluon has four components: 
[G++,G+-,G-+,G~~] = [G'',G<,G>,g'^], marked by 
the positive or the negative time branch, respectively. 
One can also express the GF in the physical representa- 
tion, where the GF is expressed by the symmetric (C), 
retarded (R) and advanced (A) components. The two 
representations are related by a unitary transformation. 
The same expressions for the SE are also valid. 

In the CTP formalism, choosing the physical represen- 
tation for the GF and SE, the Dyson- Schwinger equation 
(DSE) for G"^ reads 

V{xi)G^ixi,X2) = - J d^x'[n^{xi,x')G^ix',X2) 

+U''ixi,x')G''ix\x2)] , (2) 

G^ixi,X2)V^{x2) = - J d''x'[G'^{xi,x')n^{x\x2) 

+G^{xi,x')U^{x',X2)] , (3) 

where the differential operators T> and ^^j^e Feynman 
gauge 1^ (a = 1) are expressed by 

T^pl = 9paD';,^[A]D^;[A]+2gf'^^^F^,[A] , (4) 
V^^,^^ ^ g,,D^^^'^[A]D^^'^'^[A]+2gf'^'^'^F-,[A], (5) 

where Dj]^^[A{x)] = [dxfi +*5'^m(^)]*"' conjugate 
covariant derivative where the differential operator acts 
on the function in its left. 

In the evolution of the gluonic system one distinguishes 
different scales, which characterize quantum and soft col- 
lective motion. We introduce a mass parameter, /i, as the 
separation point of the quantum and the kinetic scale. 
In the weak coupling limit 5 <C 1, the scale of collectiv- 
ity ^ is much larger than the typical extension 
of hard quantum fluctuations ^ 1/fi. The effect of the 
classical field A on the hard quanta involves the cou- 
pling gA to the hard propagator and is of the size of 
the soft wavelength ~ 1/ (g^) . The above separation of 
scales is the basis for the gradient expansion where one 
expresses all 2-point GFs in terms of the relative y and 
the central X coordinate. Here are some typical scales: 
y^xi-x2-- X = {xi+x2)/2 - A{X) - ^ 
and F[A{X)] - gfi^. 

In order to obtain the gauge-covariant kinetic equation 
one uses the gauge-covariant Wigner function G{X, y) 
(G^, G> or G<) defined by 

G{xi,X2) = V{xuX)G{X,y)V{X,X2) , (6) 

where V{zi, Z2) — Tp exp(ig /p.^^ dz^A^) denotes a Wil- 
son link with respect to the classical background field. 
One can also define the Wilson link as a functional of 



A + Q, but it is a much more complicated case and be- 
yond the scope of this work. 

The covariant Wigner function G(X, y) transforms as 
U{X)G{X^y)U~'^{X) where only the transformation at 
a single point X is relevant. For G{xi^X2), however, the 
gauge transformation involves two points and therefore 
is not gauge-covariant. 

The DSE given by Eqs. (||) and (j^) can be expressed 
in terms of the covariant Wigner functions. To evaluate 
2?(xi)G(xi, X2) one needs to know the variation of a Wil- 
son link caused by that of its ending points. Following 
Eq. (3.15) for 6V{zi,Z2) in Ref. [note the opposite 
sign convention for g] , one finds 
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ig-SV^P^'^C! - ig-GyxFxi, 



V2 , (7) 
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-{Di,G) + {dyDxG) + {dt^G) 

3 > 1 _ 

W-:y\Fxu{dyyG) ~ ig-{dyi,G)yxFxi, 



-ig{dy^G)Ay - ig-yxFxy{Dxi^G) 
^ig^{Dx,.G)A^ + ig^G{dxi.A„) 



-ighDx.G)yxFx, - -^GA^ 
o 4 
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where we have used: Vi = V{xi,X), V2 = V{X, X2)^ G = 
G{X,y), Dxu = D^,[A{X)l = DtM,, = A^{X), 
= A^{X)A^{X) and F^, = F^,{X). The conju- 
gate expressions {ViGV2)Dl. {ViGV2)D^^^ can be 
obtained from Eqs. (|^) and (||) by taking their Hermi- 
tian conjugates and then interchanging xi and X2- Using 
these results together with Eq. (||) we derive the following 
gauge conditions for G | |26[ |: 

dxfiG^,, + ig[Gf_,^, A^] = , dyf_,G^i, = , (9) 

which are valid up to 0{gfi). The second equation is a 
transversality condition for G^,y . 

Taking the difference between Eq. (||) and Eq. (^) and 
using Eq. (||) and its conjugate expression, we derive the 
following kinetic equation for the gluon plasma: 

q-dxG^, + ig{G^,q-A^q-AG'^^ 



1 „ , 1 

C nC i?^ \ _ n (10) 



^^aivFv\{dqxG^a-^ + ^g(dq\&i-i)(li^Fi'\ 
+9{FapG^^ - G'^pFp^) = , 



where G^^ = G^^(X, q). 

The above equation is located at the collective coor- 
dinate X and is gauge-covariant under the local gauge 
transformation U{X), i.e. it transforms as U{---)U~^. 
Indeed, noting that both F^^, and G^^ are gauge- 
covariant and dq^j, does not affect /7, it is obvious that 
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the last three terms are gauge-covariant. To verify that 
the first two terms also preserve gauge covariance, we 
explicitly write down their transformation 



q ■ dxGaj 



igG^q ■ A 
-igq ■ AG'^^^ 



U{q ■ dxG%)U-^ + {q ■ dxU)G^M-' 



+ UG%q ■ 
igUG^q ■ A U-^ 
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UG^^q-dxU' 



-igUq-AG'^^U~ 
-{q-dxU)G^M^ 



With the above gauge transformations one can clearly 
see that the sum of the first two terms in Eq. (^) in- 
deed transforms asU{- ■ ■)U~^ and therefore preserves the 
gauge covariance. 

The quantum kinetic equation ( |l0|) is derived in a quite 
general and consistent manner in BG-QCD and CTP for- 
malism. No further approximations or requirements go- 
ing beyond the gradient expansion were used. We notice 
that a result similar to Eq. (|l^) was also obtained in Ref. 
p3| . There, however, the kinetic equation was derived 
by making the gradient expansion of the equation of mo- 
tion for the Wigner function (not in CTP formalism). 
In addition the derivation was made in the fundamen- 
tal (not the adjoint) color space in QCD (not BG-QCD). 
Finally the assumption was made in Ref. that the 
Wigner function is proportional to the quadratic prod- 
uct of the generators of the fundamental representation. 
The approach presented here is quite general and does 
not require any specific assumptions on the structure of 
the Wigner function. 

In the following we will compare Eq. ( pO| ) with the 
classical equation. Especially we will show the physical 
connection of the second term of Eq. (10) to the color 
precession. The classical kinetic equation for the color 
singlet distribution function f{x,p,Q) has the following 
form |,[l6[|l|,|l|, 0: 



P^.id^.-gQ^F^.dp, 
-gr'^Al{x)Q^dQ.]f{x,p,Q) 



(11) 



where is the classical color charge and a = 1, .., — 
1. The last term of Eq. ( pi] ) describes the rotation or 
precession of the color charge |||, |l^, . 

Comparing Eq. (|l^) with the quantum expression Eq. 
( p^ ) it is clear that the color singlet distribution function 
/ is replaced by the gauge-covariant Wigner function G"" 
which is a color matrix in the adjoint representation. One 
can also recognize that the first, third and fourth terms 
of Eq. (jl^) are the quantum generalization of the first 
two terms in Eq. ( ITTl) . The last term in Eq. (|o|) appears 
from the covariant operators | |27[ | and hence is not present 
in the classical equation. 

Particularly interesting is the appearance of the sec- 
ond term in Eq. (|lo[). We have seen that its presence is 
crucial to assure the gauge covariance of the Vlasov equa- 
tion. This term has an interesting physical meaning. It 
is the quantum analogue to the color charge precession 



in the classical kinetic equation. To see this more clearly, 
one can expand G'^p{X, q) with respect to the expansion 
parameter gT°'A'^^{X) from the Wilson link in Eq. (||). 
This expansion can be also understood as the result of 
the AQQ, AQQQ and AAQQ vertices. Then we have 

G%{X,q) = NoMX,q) + gT''Al{X)Ni^p.JX,q) 

+g'T''T''A;{X)AliX)N2o.p-MX, q) 
+ •••, (12) 

where T" are quantum analogues to the classical color 
charges Q"; Niap{X, q) with i = 0, 1, 2, • • • are color sin- 
glet functions. Each term of the expansion corresponds 
to an order of the color inhomogeneity in the gluonic 
medium due to its interaction with the background field. 
If the background field is refered to the soft mean field, 
its magnitude should vanish when the system approaches 
equilibrium. Then only the first singlet term survives 
m Eq. (|l|), which means the color homogeneity of the 
gluonic medium. This is somewhat similar to the multi- 
pole expansion for an electromagnetic source where the 
moments of dipole, quadrupole etc. describe increasing 
orders of spatial inhomogeneity for the electromagnetic 
charges. In weak coupling, as the lowest order approx- 
imation, we keep only the first two terms in Eq. (|l^). 
Then the second term of Eq. (noh becomes 



igiG'^^q -A-q-AGi 



-gr'^q^^A'T^dTM^^ (13) 



which reproduces the classical color precession term in 
Eq. ©. 

Since we know that Dx ~ gA{X) ~ g/z, the first two 
terms of Eq. ([lo|), i.e. q ■ dxG'^^ and the color precession 
term, are at leading order 0{g^^), while other terms are 
at subleading order 0{g^iJ?). In the vicinity of equilib- 
rium the natural scale in the system is the temperature 
T. The mean distance between particles is of the order 
of ~ 1/T, while l/{gT) characterizes the scale of collec- 
tive excitations M, 15|. For small coupling constant g 
these two scales are well separated. The covariant Wigner 
functions can be expanded around their equilibrium val- 
ues: G — + SG, where the equilibrium function G^*'-' 
is a color singlet and the fluctuation SG ~ g^G. Typical 
scales are q^T,Dx ^ g'^T, gF ~ (Dx)"^ ^ 
at leading order, only the first term of Eq. (|lO| ) survives 
and the precession term vanishes due to the color-singlet 
nature of The linearized version of Eq. (|lo| ) with 

respect to 5G corresponds to the equation formulated in 
the background Coulomb gauge in Ref. 0. 

The quantum fluctuations near equilibrium were also 
considered in Ref. in the context of the classical col- 
lisionless transport equation. It is quite natural to carry 
out the same study from our quantum approach. First, 
BG-QCD deals with the classical field and the quan- 
tum fluctuation in a systematic way. The quantum field 
plays the similar role to the field fluctuation in Ref. [ p4| . 
Second, in the quantum approach, corresponding to the 
phase-space distribution, we deal with the GFs which 
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can be expanded around their equilibrium values follow- 
ing Eq. One also needs to complete the equations 
by including the field equation D^F^^^ — (J^)- Here the 
averaged induced current (J,^) is related to the 2- and 
3-point GFs which finally depend on G. Thus it can 
be expanded according to Eq. (|lj) as well. 

The analogy and difTerences of the quantum and the 
classical Vlasov equation can also be shown by comparing 
the equations for the color moments. Corresponding to 
Eq. (|l0|) one gets 



q ■ ux 



a/3"/37 ^apFf^j) — ^ 

dxh% + gr''q-A'hl^ 



(14) 



(15) 



and 



where we defined = Tr(G^^), h^^ = Tr(r"G^^, 
/i^^ = Tr(T°T*G^^). The classical equations for color 
moments of f{x,p, Q) can be found in Ref. Compar- 
ing Eq. ( 14) and ( p^ ) with the classical expressions in Ref. 
1^, one sees that apart from the term {Faphp^ — hapFp^)^ 
which comes from the covariant operators, the quantum 



and classical equation have similar structure. The iden- 
tification of the color precession term in Eq. ( |l0[ ) is then 
st r aight forward . 

In summary, by applying the CTP and BG-QCD for- 
malism we have derived the kinetic equation for the 
gluon. The derivation is more general and transparent 
than earlier works. The kinetic equation is with respect 
to gauge-covariant Wigner function, which is a matrix 
in adjoint color space. A notable feature of our work is 
that a term is found which, as in the classical case, corre- 
sponds to the color precession of the parton. This is the 
non-Abelian analogue to the Larmor precession for par- 
ticles with the magnetic moment moving in a magnetic 
field. We see that this term is necessary to the gauge 
covariance of the kinetic equation. 
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